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Abstract 

Nagaoka's theorem on ferromagnetism in the Hubbard model with one 
electron less than half filling is generalized to the case where all possible 
nearest-neighbor Coulomb interactions (the density-density interaction V ^ 
bond-charge interaction X, exchange interaction F, and hopping of double 
occupancies F') are included. It is shown that for ferromagnetic exchange 
coupling [F > 0) ground states with maximum spin are stable already at 
finite Hubbard interaction U > Uc- For non-bipartite lattices this requires a 
hopping amplitude t < 0. For vanishing F one obtains f/c — ?■ oo as in Na- 
gaoka's theorem. This shows that the exchange interaction F is important 
for stabilizing ferromagnetism at finite U. Only in the special case X = t the 
ferromagnetic state is stable even for F = 0, provided the lattice allows the 
hole to move around loops. 
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I. INTRODUCTION 



The single-band Hubbard model was originally introduced as a model for ferromagnetism 
of itinerant electrons. ^^"^ The model is given by 



where cf^ (qv) creates (annihilates) an electron with spin a =t,i, "iT-ia = cf^c-^ is the number 
operator, and (ij) denotes nearest neighbors. This is the simplest possible correlation model 
for electrons on a lattice. However, rigorous evidence for itinerant ferromagnetism in this 
model is very limited. One of the most important results is Nagaoka's theorem, which states 
that if the Hubbard repulsion U is infinite, the ground state has maximum total spin on 
certain lattices in the case of precisely one hole. The physical mechanism behind Nagaoka's 
theorem is the following. If U = 00, the ground state of (1) is macroscopically degenerate. 
This degeneracy is lifted by the motion of the hole since it is energetically favorable for it 
to move in a background of fully aligned spins (provided the lattice allows for motion of the 
hole around loops'*). A simpler proof of Nagaoka's theorem was later given by Tasaki,^ who 
also showed that additional density-dependent interactions do not alter this result. 

Several other mechanisms leading to ferromagnetism in the Hubbard model have been 
discussed since then.^ Lieb'' proved that the ground state is ferromagnetic for bipartite 
lattices with different numbers of sites in each sublattice. Mielke® and Tasaki^ proved the 
stability of ferromagnetism for special lattices with flat bands. *° Recently, a new route 
was taken by Miiller-Hartmann** who studied ferromagnetism at low particle density in 
dimension d = 1. He included next-nearest neighbor hopping in such a way that the band 
has two minima. At low density, the on-site repulsion U generates a ferromagnetic exchange 
coupling between particles in these two pockets. 

Clearly, it is still a long way to a true understanding of itinerant ferromagnetism in 
solids. It is quite obvious that the single-band Hubbard model is not a generic model for 
ferromagnetism. So far, either the assumption of a special kind of hopping, or of [/ = 00, 
or both, were necessary to prove the stability of ferromagnetism. One may therefore ask 
if there exist other, simple mechanisms leading to itinerant ferromagnetism that are not 
contained in the Hubbard model. There are two important candidates: (i) band degeneracy, 
as it exists in 3c/-transition metals, and (ii) nearest-neighbor exchange interaction, which is 
always present in a fermionic system with Coulomb interaction. Here we discuss only the 
latter, since the effect of band degeneracy will be discussed separately.*^ 

From atomic magnetism it is known that there exists a ferromagnetic Hund's rule inter- 
action between orbitals on the same atom. A ferromagnetic spin wave function is symmetric 
under exchange of particles, and for the electron wave function to be overall antisymmet- 
ric its coordinate part must be antisymmetric itself. An antisymmetric coordinate wave 
function is zero whenever two particles are at the same position and thus minimizes the re- 
pulsive Coulomb interaction. This is the well-known mechanism leading to the flrst Hund's 
rule for atomic magnetism. In solids a similar ferromagnetic Heisenberg exchange term is 
present even between orbitals at different sites. However, since their overlap is much smaller 
than for orbitals on the same atom, this interaction may be quite small. Nonetheless, it 
will not be strictly zero. Therefore this direct exchange interaction, denoted by F below. 
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provides a natural way for stabilizing ferromagnetic states. ^"^ Of course, other features of 
the model, in particular the hopping t and the structure of the lattice, are also important 
factors concerning the stability of ferromagnetism in the ground state. 

With this in mind it is worthwhile to review the steps that originally led to the Hubbard 
model, and to retain, in a systematic way. Coulomb interaction terms beyond the on-site 
repulsion U. This is done in Sec. II, and a model Hamiltonian with all nearest neighbor 
interactions is derived. In Sec. Ill we state sufficient stability conditions for ferromagnetic 
ground states in the case of one hole in a half-filled band. In particular, it turns out that 
if the direct exchange is ferromagnetic [F > 0), and even if F = in a special case, the 
on-site repulsion U need only be larger than a finite value [/c, thereby generalizing Nagaoka's 
theorem to finite U. The details of the proof, using a method employed previously for the 
case of half-filling,^'* are deferred to the Appendix. Sec. IV contains our conclusions. 



II. DERIVATION OF THE MODEL 

Let us first review the derivation of effective models for metallic ferromagnetism. The 
general electronic model expected to describe ferromagnetic phase transitions in transition 
metals was introduced by Hubbard.* It is given by the electronic Hamiltonian*^ 



tj ijmn 

Here, cf^^ (ciaa) creates (annihilates) an electron with spin cr in a Wannier orbital a localized 
at site i. The first term describes hopping between two sites j and contains the kinetic 
energy and the ionic potential [/,o„(r). The second term describes the (screened) Coulomb 
interaction between electrons, Vee(r — r').*^ The matrix elements, expressed in the Wannier 
basis, are [h = 1) 

t" = (««| - -^V^ + U,o„{r)\ja) (3a) 
■' Am 

V"lmn = j/5|Ke(r - F ')\mfi, Tiu) (3b) 

So far no approximation was made. The Hamiltonian (2) contains infinitely many parame- 
ters. For simplicity it is therefore often assumed that the essential physics of the problem 
is captured by a single s-band, whereby all other bands are neglected. More precisely, all 
other bands are projected onto one single effective s-band. This approximation requires the 
existence of a band gap above the effective band. Then the deviation of the parameters 
tij and Vijmn from their multi-band values can be determined, in principle, by perturbation 
theory. 

The restriction to a single s-band entails considerable simplifications: orbital indices 
may be dropped in eq. (2) and (3); furthermore, all matrix elements depend only on the 
separation of the lattice sites (and not on direction). Since the matrix elements are expected 
to fall off quickly with distance, one usually retains only the first few of them. Thus hopping 
is restricted to nearest neighbor sites i and j: —t = tiy It is also natural to assume that 
U = Van is the largest matrix element of the Coulomb interaction. Keeping only t and U 
one obtains the Hubbard model, eq. (1). 
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However, there are other terms that can be of appreciable size.^ These are the two-site 
terms of the interaction: V = t^ijij, X = Umj, F = t^ijji, F' = t^ujj, where i and j are nearest 
neighbors. Keeping these terms one obtains the following single-band model: 

{ij) {ij)'^ 

+F E ^^'^^A^ + F' Ei^^^An + h-c. ) . (4) 

{ij)aa' (ij) 

Here V is the density-density interaction between nearest neighbors, X is the bond-charge 
interaction giving rise to correlated hopping, F is the exchange interaction discussed in 
the introduction (ferromagnetic in nature if F > 0), and F' represents hopping of double 
occupancies. 

While the on-site interaction U usually has the largest numerical value, the other matrix 
elements are certainly not zero. Hubbard's estimates^ for transition metals are, for example, 
U ~ lOeV, V ~ 2-3eV, X leV and F, F' ^eV, and the hopping amplitude t typically 
ranges between O.beV and l.beV. Even if nearest neighbor interactions are very small, they 
can be qualitatively important if they have different symmetries than the [/-term and thus 
can lift degeneracies. 

The model (4) was essentially derived already by Hubbard.^ Extensions of the actual 
Hubbard model (1) by some or all of the terms in eq. (4) received much attention since 
then. Eor example, Campbell, Gammel and Low^'' presented a detailed investigation of 
the phase diagram of Hnn in dimension o? = 1, and discussed the relative magnitude of 
its parameters for real materials. On a mean-field level, the effect of the terms in (4) on 
the stability of ferromagnetism was studied by Hirsch.^® Eurthermore, exact solutions are 
possible in the special case of X = t. In this case the number of doubly occupied sites is 
a conserved quantity, and the exact ground state solution can be obtained in a wide range 
of parameters. ^^^^^ Eor X = t and V = F = F' = the model was recently solved exactly 
in one dimension, ^"^'^'^ while for X = t = —V = F = F' a solvable supersymmetric model is 
obtained. ^^'^^ The case X = t will play a special role in our analysis, too. 

Criteria for the stability of ferromagnetic ground states of the Hamiltonian (4) were 
recently derived for the case of half-filling (one electron per site).^'*'^^^^^ The ferromagnetic 
states are then found to be insulating. To gain insight into the more general problem of 
itinerant ferromagnetism we will now investigate a half-filled band with one hole, as in 
Nagaoka's work."* Thus we consider a finite lattice with L sites and fix the total number of 
particles at N = L — 1. The number of nearest neighbors is denoted by Z. We consider 
lattices with at least Z nearest neighbor bonds between any subset of lattice sites and the 
set of remaining sites. Eor example, all crystal lattices with periodic boundary conditions 
fulfill this requirement. 

III. FERROMAGNETIC GROUND STATES 

The Hamiltonian Hnn commutes with the total spin S = J2i Si, where = 
I cf^'^'aa'^ia' ^^^^ Pauli matriccs. The eigenvalues of are denoted by S'(5'-|-l). 

In the following we will be concerned only with saturated ferromagnetic states with largest 
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possible eigenvalue Smax = N/2 = [L — l)/2. There are 2Smax + 1 = iv such states with the 
same energy eigenvalue. 

We are interested in the following question: Under which circumstances do the ground 
states of Hnn have maximum spin? For the pure Hubbard model ii. V = X = F = 
F' = 0), Nagaoka's theorem"* states that for U = oo,t<0(t7^0if the lattice is bipartite) 
the ground states have S = Smax- This statement can be generalized to arbitrary density- 
density interaction V ."^ These results require the lattice to have "loops" as discussed in 
the introduction. For example, this is the case for the square, triangular, simple cubic, 
body-centered cubic, face-centered cubic, and hexagonal close packed lattice, but not in one 
dimension or on the Bethe lattice.'* 

The main result of this paper is a generalization of Nagaoka's theorem to finite values 
of the Hubbard interaction U . In the Appendix we derive the following conditions for 
ferromagnetic ground states: 

The ground states of H^n with one hole (i. e. N = L — 1) have maximum total spin 
S = Smax = in the following cases: 

Case 1: On any lattice, if F > 0, t <,0 and (a) X ^ t and U > Uj:^\ or (h) X = t and 

Case 2: On lattices with loops, if X = t < 0, F = 0, and U > [/^^ . 
In both cases t > is allowed if the lattice is bipartite. 

These results are summarized in Table L The constants U^^'^ and [/^^ are given by 



up = Z{2\t\ + 
U^ = Z(2\t\ + 



V - F -2\t\ 
F 

V 2t 

2 ' ' 



+ 



F 



+ 



F' 



(X-t) 



F 



+ \F' 



(5a) 
(5b) 



Hence, if F > 0, ferromagnetic ground states are stable on any lattice for U larger than 
a finite critical value. For F — > 0"*" we have U^'^ — > cx3, thus yielding Nagaoka's condition for 
the pure Hubbard model. This shows that the Heisenberg interaction F, which is neglected 
in the Hubbard model, provides an obvious mechanism for stabilizing ferromagnetic ground 
states at finite U. Note that since X and t are expected to be of the same order of magnitude, 
the sensitive dependence on F, due to the term may cancel from Uj:^\ and values of 

the order of Uc ~ 12eV are possible. The dependence of Uc on t, F is depicted in Fig. 1. 
The case X = t is special, since in this case the stability of ferromagnetism can be achieved 
either by F > 0, or by F > and t < if the lattice has loops. 



Case 




Condition on U 


Condition on lattice 


Condition on t 


la 


F > 0, X / i 


u > 


any lattice 


bipartite lattice: 
non-bipartite lattice: 


t arbitrary 
t < 


lb 


F >0, X = t 


u > uP 


any lattice 


bipartite lattice: 
non-bipartite lattice: 


t arbitrary 
t < 


2 


F = 0, X = t 


u > uP 


lattice with loops 


bipartite lattice: 
non-bipartite lattice: 


t / 
t < 



TABLE I. Sufficient conditions for ferromagnetic ground states with one hole. 
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FIG. 1. Critical value Uc vs. exchange interaction F for different t, V, X , and F' = 0. For 
U > Uc the ground state is ferromagnetic. See Table I for details. 

The critical couplings Uj:^'^ and [/^^ are sums of terms, each of which corresponds to 
a typical energy scale. This means that the on-site interaction U has to be larger than 
the energy describing the paramagnetic state (bandwidth ~ Z\t\), as well as the threshold 
energies for the onset of a charge-density wave or phase separation (~ Z\V\), ry-pairing 
superconductivity^^ (~ and a spin-density wave (~ [X — ty / F). Note, however, that 

these terms do not enter separately, but appear in combinations, i. e. the effects interfere as 
should be expected. 

We remark that the above conditions are sufficient conditions. The occurrence of ground 
states with maximum spin outside the above parameter region is not ruled out. For example, 
it may still be possible to find ferromagnetic ground states even for X ^ F = and 
Uc < oo. It is interesting to note that the values for Uc appearing in eq. (5) are the same 
as those for the case of half-filling (no hole).^"* The bound Uc for half-filling was recently 
improved by Schadschneider and de Boer.^^ 

As shown in the Appendix, the ferromagnetic ground states are the same as those dis- 
cussed by Nagaoka,"* i. e. the wave function with = Smax corresponds to a band filled 
with spin-up electrons, with the hole at the top of the band, as illustrated in Fig. 2. If t < 
the band maximum is at the origin, and the corresponding wave function is I^/jq) = aQ||"t^), 
where a^^^ = -^J2^exp{i\il{.i)c-^, and I'D = Ylicf^\0) = Ilk^ktlO) the filled spin-up band. 
For t > a bipartite lattice is required and the band maximum is at wave vector Q defined 
by eic+Q = — ^k, e. g. for a hypercubic lattice we have Q = (tt, tt, . . . tt). In this case the wave 
function is Itpo) = cIqJ^). All 2Smax + 1 ground states can be obtained from Itpo) by global 
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FIG. 2. The figure illustrates (for d = 1) that the ferromagnetic ground state corresponds to 
a band filled with spin-up electrons (indicated by arrows), from which a single electron has been 
removed at the top of the band (open square), (a) i > (with a hole at A; = tt); (b) i < (with a 
hole at A; = 0). 

SU{2) rotations, i. e. the state 

\^pM) = {S^-iSy)^\^po) (6) 

has = Smax — where M = 0...iv — 1. The states {|V'm)} have the same energy 
eigenvalue. 

If F > 0, our criteria permit ferromagnetism in any spatial dimension. In particular they 
hold in dimension d = 1. This is not in conflict with the Lieb-Mattis theorem^'' which rules 
out ferromagnetism for any continuum model with A;^-dispersion and symmetric interaction 
potential in o? = 1. In a Bloch basis a A;^-dispersion corresponds to inflnitely many bands 
as in (2), whereas the Hamiltonian H^n in (4) is obtained by a projection onto one of 
these bands, as discussed in Sec. II. Hence the model (4) is quite different from the initial 
multi-band Hamiltonian H in (2), in particular there does not exist a potential energy in a 
continuum model that gives rise to exactly the terms appearing in the truncated Hamiltonian 
Hnn- Therefore the Lieb-Mattis theorem does not apply to the single-band Hamiltonian (4) 
in 0? = 1. In the Appendix of Ref. 27 Lieb and Mattis proved a second theorem, using 
the occupation number formalism, which precludes ferromagnetism in o? = 1 in the case 
of next-neighbor hopping and purely density-dependent (i. e. site-diagonal) interactions. 
Therefore this second theorem also does not apply to the Hamiltonian H^n^ eq. (4), where 
the non-diagonal exchange interaction F > (which is always present as part of the Coulomb 
interaction) plays a crucial role. The role of non-diagonal interactions for the stability of 
ferromagnetism was already discussed by these authors.^'' 

IV. CONCLUSIONS 

We presented a generalization of Nagaoka's theorem to a Hubbard model with all nearest- 
neighbor interactions. For this model, with one hole in a half-flUed band, we derived rigorous, 
sufficient conditions for the stability of saturated ferromagnetism in the ground state. The 
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ferromagnetic ground state is found to be stable on any lattice provided the next-neighbor 
exchange F is ferromagnetic and the Hubbard repulsion U is larger than a critical value 
[/c < CX3, with Uc ~ 1/F for F — > 0"*". If F = 0, only the special case X = t (unlikely to be 
fulfilled exactly in real materials) yields a ferromagnetic ground state on lattices with loops. 

The ferromagnetic ground state is an itinerant state with non-zero kinetic energy. The 
proof of its stability cannot be easily extended to doping beyond a single hole. Of course, 
this would be highly desirable since a single hole is irrelevant in the thermodynamic limit. 
However, the ground states of the model with next neighbor interactions are difficult to 
obtain for finite hole densities, since in this case simple eigenstates of the Hamiltonian are 
not known, such that the present methods can be applied only in special cases. 

To be able to explain ferromagnetism in more detail, it is clear that more ingredients 
are needed than those contained in the single-band Hubbard model with nearest neighbor 
hopping. Of greatest interest is the case of band degeneracy where the present methods can 
be applied, too.^^ 

The authors would like to acknowledge valuable discussions with D. Baeriswyl, E. H. Lieb, 
D. C. Mattis, E. Miiller-Hartmann, and A. Schadschneider. M.K. and R.S. gratefully ac- 
knowledge scholarships of the Studienstiftung des Deutschen Volkes. This work was sup- 
ported in part by the Sonderforschungsbereich 341 of the Deutsche Eorschungsgemeinschaft. 

APPENDIX: DETAILS OF THE DERIVATION 

To derive the sufficient conditions for ferromagnetic ground states stated in Sec. HI we 
rearrange the Hamiltonian (4) as a sum of positive semi-definite terms. This enables us to 
construct a lower bound on the ground state energy. If this lower bound coincides with the 
eigenvalue of a trial state, this state is an exact ground state of (4).^®'^^ Here we consider a 
wave function of the form 

IV'o) = ^ajCjfl t) , 

i 

where | 'f) = JJ- cf^\0) is the state with all sites occupied by spin-up electrons. The coefficients 
ttj ^ will be determined below. Itpo) is a state with N = L — 1 particles. Since it is an 
eigenstate of with eigenvalue Smax^ it is a representative of the subspace with S = Smax- 
It is sufficient to consider a state of this form, because whenever this state is a ground state, 
so are its (25* -\- 1) = L global SU{2) rotations. 
We introduce the operators 

P^Ja = (1 - rl^-a)ic,^ + AlC^-^)(l - %_^) 

Here Ai = sgn(t), A2 = sgn(X — t), A3 = sgn(F' — \X — and a 7^ is an arbitrary, 

real parameter. Eurthermore we introduce 
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_ I e«ej + 3d,dj + p,dj + pjd, if V > 2\t\ + {F + a^\X - t\)/2 
~\^{p.-PjY + 2{le, + d,e,) ifV <2\t\ + {F + a^\X -t\)/2 ' 

where = (1 — ?^it)(l — J^ij,), Pi = (j^if — ^^4)^, di = nifn^ are the projectors onto an empty, 
singly and doubly occupied site, respectively. It is straightforward to verify that H^n can 
be written as 



Hnn = e [1^1 E(4.% + QtQ.,.) + \x-t\ i+4 + 

+1 v^l E + f'^tf^^i - E s.s, , (Ai) 



where 



F' 


= F' 


F 


= F 


V 


= V 


U 


= U 



\x -t\ 



a2 



(A2a) 

a'\X-t\, (A2b) 

F + a^\X-t\ , , , 
^ --m. (A2c) 

(A2d) 



a2 



and an overall constant was dropped. 

Let us consider the terms in (Al) one by one. The first term is positive semi-definite. 
Since A-ltpo) = 0, B-ltpo) = 0, and Q-^ltpo) = 0, the action of this term on Itpo) is given by 

(uV {ij)ma 

= 1^1 E(«»' + -^i«j)(c,t + -^is-t)lt) • 

Thus, if tti = — sgn(t)aj for all nearest neighbors we find that Itpo) has zero eigenvalue 
and hence is a ground state of this term. Let us consider the case t < first. Then we must 
have tti = ttj, and the resulting wave function is 

iV'o) = 4f E^.fl t) = «k=otl t) , 



where Oj^^ = J2 - exp{i]iRi)c-^. Note that the filled spin-up band can be expressed as | "j^) = 

rik ^kfl^)- The (non-interacting) band structure is given by ek = — ^l^(ij) exp(zk(Rj- — Rj)). 
Thus for t < the hole at A; = is created at the band maximum^ where ek=o = Z\t\. This 
is shown in Fig. 2a. 

Next consider t > 0. Then = —a^ must hold for all nearest neighbors (ij). This is only 
possible if the lattice is bipartite, i. e. if any two nearest neighbors are located on different 
sublattices A and B. In this case the wave function becomes 
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IV'o) = ^(E^.t - E^.t)lt) = «k=Qtlt) • (A3) 

Here Q is defined by exp[i(^(R- — Rj)) = —1 for nearest neiglibors or equivalently by 
ek+Q = — Ck- Tliis implies that the band maximum is at k = Q for t > 0, and this is the 
state in which the hole is created, as shown in Fig. 2b. 

From now on we choose the coefficients as just described, Itpo) thus being a ground 
state of the first term in (Al). If the lattice is not bipartite, this requires t < 0. (The trivial 
case t = can also be included here.) 

Turning to the positive semi-definite 0-term in (Al), we observe that for V > it 
annihilates |?/'o), which in this case is a ground state of this term. Now consider V < 0, in 
which case Itpo) is an eigenstate of Y^^ij^ 0- with eigenvalue Z/2. We now show that this 
is the lowest possible eigenvalue for N = L — 1. Consider states with a fixed configuration 
{pi = 0,1} of P singly occupied sites, P = J2iPi- The possible (integer) values of P are 
0< P < L. 

(i) P = L. This case is impossible for N = L — 1. 

(ii) 1 < P < L — 1. In this case we have 

{ij) {ij) {ij) 

i. e. 20 is bound from below by the number of bonds (ij) with pi ^ py Consider the set of 
sites with pi = 1. By our definition of a lattice at the end of Sec. II, these sites are connected 
to at least Z sites with pj = 0. Hence > Z/2. 

(iii) P = 0. Then, since N = P — 1, the number of lattice sites L must necessarily be odd, 
and there must be empty and doubly occupied sites. (There are no singly occupied 
sites.) Consider states with a fixed configuration of doubly occupied sites {di = 0, 1}. We 
have 

> 2^((/,(l - d,) + d,{l - di)) = 2E(1 - ^d.d,) . 

Fxcept for a factor of 4, this is the same problem as above, with {pi} replaced by {di}. 
Therefore in this case > 2Z > Z/2. 

Thus, summarizing the cases (i)-(iii), we obtain > Z/2. Hence Itpo) is a ground state 
of the 0-term for any since it is always an eigenstate with the lowest possible eigenvalue. 

Finally, Itpo) is clearly a ground state of the remaining terms in (Al) if [/ > and F > 0, 
since it has no doubly occupied sites and maximum spin. 

So far we proved that Itpo) and its global SU{2) rotations are among the ground states 
of Hnn if U > and F > (and t < if the lattice is not bipartite). Clearly, if F > 
(and U > 0) these are the only ground states, since then only states with maximum spin 
minimize the F-term. This will be used below, when we prove that all ground states of H^n 
have S = Smax in the cases listed in Table I. 

Case la: X ^ t, F > 0, U > U^/l In this case we choose = [1 - e)F/\X - t| > 
with < e < I to be specified later. Then U in eq. (A2) becomes 
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U = U - Z{2\t\ + \V - F - 2\t\ + ^ 
Using (1 — e)~^ > 1 + 2e we obtain the bound 



2(X - t 



i2 



+ max(F, -F')). 



U > U - Z{2\t\ + 



V - F -2\t\ 



+ max(F', ^ ^ ^ - F') + e(— + 



F ' ' '2 F 



= [/-[/W-||(F^ + 8(X-tf), (A4) 

where Uj:^'^ is defined in eq. (5). Since U — Uj:^'^ > by assumption, it can be seen from 
eq. (A4) that U is positive if e is chosen small enough. For example, 



e = mm 



F{u-u^'y) 1 

Z(F2 + 8(X-t)2)'2 



indeed yields U > 0. Furthermore, from eq. (A2), F = eF > 0. Hence the present choice of 
yields U > and F > 0. Therefore only states with S = Smax are ground states of H^n- 
Case lb: X = t, F > 0, U > U^, with as defined in eq. (5). In this case 

U = U — [/p^ > and F = F > 0. Again, only states with S = Smax are ground states of 
Hnn- 

Case 2:X = t<0^F = 0^U > Uj:'^\ for a lattice with loops. (In the case of a bipartite 
lattice, a phase transformation on every other site may be applied to include the case t > 0.) 
We know that Itpo) is a ground state of Hnn for U > and that it is lower in energy than any 
state with doubly occupied sites. Hence no ground state of Hnn can have doubly occupied 
sites. Among the states without double occupancies those with S = Smax are the lowest 
in energy; this follows from the P-term in (Al) by a proof completely analogous to that of 
Tasaki.^ In Nagaoka's basis, 

|z,((7l...(7,_i(7, + i •••(Jl)) = (-1)'C+^^ • • -ctl.a.-iCHl.^^ + l •••ctaJO), 

all off-diagonal matrix elements of the P-term are negative. Furthermore, it is represented 
by an irreducible matrix since it is assumed that the lattice allows the hole to move around 
loops. It follows directly from the Perron- Frobenius theorem"^° that, in every sector with 
fixed S*^, the ground state is unique and is given by a linear combination with strictly positive 
coefficients of the basis vectors. These ground states are just the states {|V'm)} of eq. (6) 
and all have S = Sm ar- 
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